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ABSTRACT 



■ We have investigated the linear polarisation of Lya transferred in a static, 

m 

dustless, and thick plane-parallel slab that is uniformly filled with neutral hy- 



drogen. Because the scattering phase function is characterised by the dipole 



. type for wing scatterings of Lya, which is in high contrast with more isotropic 

■ scattering phase function associated with resonant scatterings, the polarisa- 

tion of the emergent Lya is dominantly controlled by the number of wing 

(N ! 

, scatterings before escape from the slab. When the scattering medium is ex- 

. tremely thick, photons emergent in the grazing direction are polarised in the 

I ■ direction perpendicular to the slab normal up to 12%, which coincides with 

o 

!h ' the limit attained by the Thomson scattered radiation in a semi-infinite thick 

\ medium investigated by Chandrasekhar in 1960. The number of wing scat- 

terings before escape is sensitively dependent on the product of the Voigt 

■ parameter a and the line-centre optical depth To, with which the variation of 

the degree of polarisation can also be explained in comparison with the Thom- 
son scattering case. The polarisation direction is parallel to the slab normal 
when aTQ < 10 3 , while it is perpendicular when aro > 10 3 . Our simulated 
spectropolarimetry of Lya from a highly thick medium is characterised by 
negligible polarisation around the line centre, parallel polarisation in the near 
wing parts and perpendicular polarisation in the far wing parts, where the 
polarisation direction is measured with respect to the slab normal. The polar- 
isation flip naturally reflects the diffusive nature of the Lya transfer process 
in both real space and frequency space. We also discuss the beaming effect 
that is associated with the parallel polarisation in the near wing parts of the 
emission. 

Key words: line: profile - radiative transfer - polarisation 
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2 Lee and Ahn 

1 INTRODUCTION 

In the ultraviolet band there are many resonance lines arising from the S1/2 — > Pi/2,3/2 
transitions including CIVAA1548, 1551, NVAA1238, 1241, MgIIAA2798, 2800, and LyaA1216. 
These lines are found to be prominent and conspicuous in many line-emitting objects in- 
cluding quasars, starburst galaxies, and planetary nebulae because they are major coolants 
in plasmas with temperature T ~ 10 4 — 10 5 K. The atomic physics associated with the 
S1/2 - ► Pi/2,3/2 transitions provides a nice and simple example of the quantum interferences 
between the excited levels. 

It is well-known that the resonance scattering between two levels with J = 1/2 yields 
isotropic outgoing radiation that is completely unpolarised. This is in high contrast with the 
fact that a 90° resonance scattering of an unpolarised photon associated with a J = 1/2 — > 
J = 3/2 transition yields the maximum degree of polarisation 3/7 (e.g. Stenflo 1980; Lee, 
Blandford & Western 1994). Stenflo (1976a, 1976b, 1980, 1994, 1996) and Lee & Blandford 
(1994) studied the polarisation of resonance lines. 

Interesting phenomena occur when the incident radiation is scattered in off-resonance 
regimes, where both the polarisation and the angular distribution of the outgoing radiation 
are sensitively dependent on the wavelength of the incident radiation, as is illustrated by 
Stenflo (1980). One general result is that scattering in far off- resonance is characterised by 
the classical Rayleigh scattering phase function. 

Being the most prominent emission features with a very small fine-structure level-splitting 
in the excited states of hydrogen, Lja should be treated with more emphasis and care. Due 
to the small fine-structure level-splitting of hydrogen, wing scatterings of Lja are charac- 
terised by the classical Rayleigh scattering phase function, which is the same phase function 
associated with the Thomson scattering. Hence, when wing scatterings of Lya become im- 
portant in a thick medium of hydrogen, we naturally expect that the Lja line transfer may 
be described in a way similar to the transfer of continuum photons in a Thomson scattering 
medium. 

Chandrasekhar (1960) investigated the polarisation of photons transferred in a very 
thick electron gas, where he showed the degree of polarisation reaches 11.7% for photons 
emergent in the grazing direction of very thick and plane-parallel media. Phillips & Meszaros 
(1986) also investigated in a numerical way the radiative transfer in electron clouds, for both 
optically thin and thick cases. They found that the photons are beamed to the direction 
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normal to the slab, and obtained 11.7% of polarisation confirming the work of Chandrasekhar 
(1960). 

The vertical column density A/hi of normal galaxies and dwarf galaxies ranges A^hi = 
10 19 — 10 22 cm~ 2 . The line-centre optical depth is related with the HI column density N m by 
— 1/11 rr-i/2 \ AT HI " 

where T 4 is the temperature of the medium in units of 10 4 K. This implies that the Lya line- 
centre optical depth r = 10 6 — 10 9 for most galaxies, and therefore one needs to deal with an 
extremely thick medium. It appears that many galaxies found in the early universe exhibit 
asymmetric Lya emission, which is supposed to be emergent from media with these high 
optical depths. Therefore, if the scattering medium possesses any geometrical or kinematical 
anisotropy, we may expect that the Lya emission from these objects will be polarised. 

A preliminary investigation of the polarisation of Lya emission in a thick medium has 
been presented in Lee & Ahn (1998). They considered an anisotropically expanding and 
optically thick neutral medium with a Hubble-type flow, and found that the emergent Lya 
photons are linearly polarised with a degree of polarisation up to 10% in the direction 
parallel to the slab normal. They also applied their calculation to a hemispherical shell 
partly obscured by an opaque component such as a disc, and showed that about 5% of 
polarisation may develop. 

In this paper, we calculate the polarisation of Lya in more detail with the particular 
attention on the polarisation direction as a function of the optical depth. Section 2 describes 
the configuration of our model and the Monte Carlo code. In section 3, we show the re- 
sults, and explain microscopic processes of scatterings that lead to the development of Lya 
polarisation, which is followed by the summary section. 



2 MODEL 

In a very thick medium, wing scattering plays a very important role in the transfer of res- 
onance line photons, which is similar to a diffusion process. Previous researchers (Avery & 
House 1968; Adams 1972; Harrington 1973; Neufeld 1990) introduced the diffusion approx- 
imation where only wing scatterings take place during the transfer. This approximation is 
very efficient for the extremely thick cases with ar > 10 3 , where a is the Voigt parameter 
defined as a = T/AttAud, the ratio of the natural width T and the Doppler frequency width 
Au D = vo(vth/ c ) with v t h and u being the thermal speed and the line centre frequency, 
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4 Lee and Ahn 

respectively. We introduce the dimensionless parameter x that measures the frequency shift 
from the line centre in units of Az/p, 

x = (v-vq)/Av d . (2) 

However, Lya photons experience both core scatterings and wing scatterings and may 
alternate between these two types of scattering during the line transfer. A large number of 
core scatterings change the polarisation vector in a random way yielding an isotropic and 
unpolarised radiation field. Therefore, we have to consider not only wing scatterings but 
also core scatterings in order to compute accurately the polarisation of Lja emergent from 
a very thick H I medium. 

Ahn, Lee, & Lee (2000, 2001, 2002) developed a Monte Carlo code which can deal with 
the Lyct line transfer in thick media. They first dealt with moderate optical depths, because 
the diffusion approximation is not valid in this optical depth regime. They developed an 
accelerating scheme for the code in order to improve the computing speed, which enabled 
one to calculate the line transfer in extremely thick media. The code also accurately treats 
the partial frequency redistribution and is equipped with a subroutine that differentiates 
between resonant core scatterings and non-resonant wing scatterings, providing different 
scattering phase functions according to the type of scatterings. 

In this work, we use the same Monte Carlo code, of which the description can be found in 
our previous work. Our code treats the scattering of Lja photons in a manner very faithful 
to the atomic physics associated with the fine-structure of hydrogen. The level-splitting of 
2Pi 3, the excited state of Lja transition, is 10GHz, which amounts to the Doppler width 
of 1.34 km s" 1 relative to Lja. Because this is much smaller than the thermal speed of a 
medium with T > 100 K, one normally neglects the 2Pi | level-splitting, treating it as a 
single level. Even when T < 100 K, the level-splitting effects may be neglected in the cases 
where the line centre optical depth of the scattering medium is very large. This is because the 
escape of line photons can only be made in considerably far wings where the level-splitting 
is negligible. 

Since we are particularly interested in the development of polarisation of Lya, we provide 
the detailed description of the code about the angular redistribution and the polarisation 
of the scattered radiation. We adopt the density matrix formalism to describe the angular 
distribution and polarisation of the scattered Lja, where the density operator is represented 
by a 2 x 2 Hermitian matrix p. The usual relation with the conventional Stokes parameters 
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I,Q,U and V is given by 

I+Q I-Q 
Pn - — 2 — ' P22 ~ — 2 — 

P12 = — ^ — ' p21 = Pl2 ^ 
Fig. 1 shows the model configuration adopted in this paper. We consider a static, uniform, 
and optically thick HI slab with the total optical thickness 2r along the z— axis that is chosen 
to be the slab normal. The completely unpolarised Lya source is located at the origin of 
the coordinate system and p is defined as the cosine of the angle between the wavevector 
of the emergent photon and the slab normal. We will not consider the circular polarisation 
because of the azymuthal symmetry assumed in this work. Therefore V = and p may be 
regarded as a real matrix. 

In our Monte Carlo code, we only consider three types of scattering of Lya, which are 
resonance scattering with S1/2 — > P1/2, resonance scattering with S1/2 — > P3/2 and wing 
scattering far from both resonance transitions. At each scattering event, we compute the 
probability which type of scattering has occurred. The procedure for the computation of 
this probability is described in detail in our previous work (Ahn et al. 2000). 

The incident radiation is described by the density matrix p and the wavevector k = 
(sin 9 cos 0, sin 9 sin 0, cos 9), and the scattered radiation is described by the corresponding 
primed quantities p' and 

k' = (sin 9' cos 0', sin 9' sin 0', cos 9'). 

We choose the polarisation vectors 

ei = (— sin0, cos0, 0) 

e 2 = (cos 9 cos 0, cos 9 sin 0, — sin 9) (4) 

for an incident photon and correspondingly e'-^e^ for the scattered photon. 

With this choice of the polarisation basis and the scattering geometry shown in Fig. 1, the 
polarisation will develop either in the direction parallel to the slab normal (the e 2 component) 
or in the direction perpendicular to the slab normal (the eicomponent). Therefore, we may 
expect that the polarisation of the emergent radiation will be described completely by P = 
Pu — P22, the difference of the diagonal elements of the density matrix, which we will call 
the degree of polarisation. Hence in this convention, the positive degree of polarisation 
means polarisation in the direction perpendicular to the slab normal, and the negative 
polarisation is polarisation in the direction parallel to the slab normal. In this work, we refer 
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6 Lee and Ahn 

the polarisation direction with respect to the slab normal vector which is the symmetry axis 
of the scattering medium. 

When a wing scattering occurs, the density matrix associated with the scattered radiation 
is given by 

Pn = Pn cos 2 A0) — p 12 cos0 sin2A0) 

+P22 sin 2 A0cos 2 
p' 12 = -pn cos0'sin2A0 + pi 2 (cos0cos0'cos2A0 

+ sin 9 sin 9' cos A0 (5) 

+P22 cos 9(— sin 9 sin 9' sin A0 

— - cos cos 0' sin 2A0) 

P22 = Pn cos2 @' s i n2 ^</> 

+P12 cos 9' (2 sin sin 9' sin A0 

+ cos 9 cos 0' sin 2 A0) 

+P22(cos cos 6 1 ' cos A0 + sin 6* sin 9') 2 , 
where A0 = 0' — 0. 

When a given scattering is a resonance scattering associated with the transition Si/ 2 — > 
P3/2, then the density matrix is obtained by the following relation 

Pn = Pn(5 + 3 cos 2A0) 

+P22 [((5 - 3 cos 2A0) cos 9 2 + 2 sin 9 2 ] 

— P126 cos 9 sin 2 A0 
P12 = 3pn sin 2 A(f> cos 0' 

+6pi2 (cos 9 cos 0' cos 2 A0 

+ sin 9 sin 0' cos A0) (6) 
+3p22 cos 0(— 2 sin sin 9' sin A0 

— cos 9 cos 0' sin 2 A0) 

p' 22 = p 11 [(5-3cos2A0)cos 2 0' + 2sin 2 / ] 
+P22 [(5 + 3 cos 2 A0) cos 2 9 cos 2 9' 
+2 cos 2 sin 2 9' + 12 cos A0 cos cos 0' sin 9 sin 0' 
+2 cos 2 9' sin 2 + 8 sin 2 sin 2 0')] 
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Lya Polarisation in Thick H I Media 7 




Figure 1. Configuration of the model adopted in this work. We consider the static and uniform slab with the Lya source 
at the origin denoted in the figure by O. The vertical thickness of the slab is 2-ro, and we define fi as the cosine of the angle 
between the wave vector of the emergent photon and the slab normal. 

+P12 (6 sin 2 A0 cos 2 9' cos 9 
+2 sin A0 cos 9' sin 9 sin 9') . 

Finally, when the scattering is resonant with the transition S1/2 — > P1/2, then the density 
matrix for the scattered radiation is set to be 

P'u = P22 = \, P12 = P21 = 0, (7) 

which corresponds to the istropic and perfectly unpolarised radiation. 

The angular distribution of the scattered radiation for an incident photon with k and 
p is given by the trace of p' (9 1 , <fr' ', 9 , 0) , as is shown by Lee et al. (1997). Therefore, once 
the wavevector k of a scattered photon is chosen from p' in accordance with the scattering 
type, the polarisation is also determined at the same time. It is also noted that the degree 
of polarisation is computed only after the new density matrix is normalised so that it has a 
unit trace. 



3 RESULTS 
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Figure 2. Polarisation of the emergent Lye? photons originated from the unpolarised source located at the midplane of the 
slab illustrated in Fig. 1. We show five cases with different a to. The solid line in each box represents the limiting behaviour for 
the Thomson scattering. 

3.1 Degree of polarisation 

We calculate the degree of linear polarisation of Lya photons emerging from an optically 
thick slab with unpolarised sources distributed in its central plane. In Fig. 2, we show the 
linear degree of polarisation for various optical depths according to /x that represents the 
observer's line of sight with respect to the slab normal. In the calculations, we consider the 
cases where optical depths to the normal direction of the slab are tq = 2 x 10 4 , 5 x 10 4 , 2 x 10 5 , 
5 x 10 5 , and 2 x 10 6 . Since we fixed a = 1.49 x 1(T 2 , ar = 3 x 10 2 ,8 x 10 2 ,3 x 10 3 , and 
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8 x 10 3 , respectively. In the figure, the dots with la error bars represent our results and the 
solid lines represent the result for continuum photons transferred in a Thomson scattering 
semi-infinite electron cloud (Chandrasekhar 1960). We can see in the figure that the case 
with ar > 10 3 shows good agreement with the limiting case of the Thomson scattering. We 
also see that the polarisation curves gradually converge to that of the Thomson scattering 
with a very large optical thickness as ar increases. 

When aro > 10 3 , the Lya transfer processes are dominated by wing scatterings accom- 
panied by the spatial transfer and Lya photons experience a number of successive wing 
scatterings just before escape (Ahn et al. 2002). The number of last successive wing scatter- 
ings is determined by the scattering optical depth, t w at the characteristic frequency of the 
emergent photons, given approximately by 

1~w /=%Sl (8) 

where x s = (aro) 1 / 3 . They also showed that these are valid only for ar > 10 3 , when the 
diffusion approximation can be safely applied. 

For the cases considered in this section, the characteristic wing optical depths are t w = 
3.8,5.1,8.1,11.0,17.5, respectively. According to Phillips & Meszaros (1986), the limiting 
behaviour is almost achieved when r e > 10. The polarisation of Lya exhibits the same 
limiting behaviour to that of the Thomson scattering when ar > 10 3 or t w > 10. 

3.2 Spectropolarimetry 

We are interested in spectropolarimetric observations for many astronomical objects. In 
Figs. 3-7 we show our results of the accelerated Monte Carlo calculations for the polarisation 
of Lya emerging from an extremely thick hydrogen slab. The total number of photons in 
each simulation is 1.6 x 10 6 , and the frequency bin is set to be Ax = 1. For r = 2 x 10 6 
we set Ax = 2 as an exception. All the quantities shown in the figures have been obtained 
after averaging over all emergent angles 9 = cos -1 //. 

In the top panel of Figs. 3-7, we show the emergent fluxes, in the middle panel the degree 
of polarisation Pf, and in the bottom panel the characteristic optical depth, t w . Treating the 
radiative transfer in a Thomson scattering medium as a random walk process, we convert N w , 
the average number of successive wing scatterings just before escape, to the characteristic 
optical depth t w = y/N w . 

It is immediately seen that negligible polarisation develops near the line-centre. This is 
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Figure 3. Simulated spcctropolarimetry of Lya emission for the unpolariscd source at the central plane of a slab with the 
vertical optical depth from the centre of the plane to = 2 X 10 4 and the Voigt parameter a = 1.49 X 10 — 2 . The top panel shows 
the emergent flux, the middle panel shows the degree of polarisation, and the bottom panel shows the effective wing optical 
depth = \/N w where N w is the number of successive wing scatterings just before escape. All the quantities shown in the 
figure have been obtained after averaging over all emergent angles = cos -1 ft. 
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Figure 4. Simulated spcctroplarimetry of Lya from a slab with ro = 5 X 10 4 . All other quantities represent the same as in 
Fig. 3. 
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Figure 7. Simulated spcctroplarimetry of Lya from a slab with to = 2 X 10 6 . All other quantities represent the same as in 
Fig. 3. 

because line photons with frequencies near the line-centre escape through a single longest 
flight after a large number of core scatterings near the surface, where the local core scatterings 
isotropise the radiation field. The core-wing boundary frequency is approximately given by 
Eq.(8) with t w — 1, or x — ypH. Hence our results in the figure show that Pf — at 
< x < v^F. 

The second point is that the photons filling the near wing part of Lya emission are 
positively polarised, which means the electric vectors of those escaping photons tend to lie 
in the direction perpendicular to the slab normal. The degree of polarisation in this part is 
about 1%, which is smaller than the maximum value 11.7%, because they are averaged over 
emergent angles. It is very notable that the degree of polarisation becomes negative at far 
wing frequencies, which means that the polarisation direction flips from perpendicular to 
parallel as the frequency shift Ax from the line-centre increases. 

In order for wing photons to acquire scattering numbers larger than N w = r^, they 
spatially diffuse much farther than r in the stage of their final series of wing scatterings. 
According to Ahn et al. (2002), a line photon smears out towards the slab surface via wing 
scatterings, and therefore the photon diffuses into the grazing direction of the slab acquiring 
a sufficiently large frequency shift. This means that the scattering medium is effectively 
very thin in the normal direction from the view point of far wing scattering Lya photons. 
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Lya Polarisation in Thick H I Media 13 

Therefore, the scattering plane associated with the escaping photons nearly coincides with 
the slab plane, leading to the development of polarisation in the direction parallel to the slab 
normal. This explains the negative polarisation in the far wing part illustrated in Figs. 3-7. 

On the other hand, the near wing part of the emergent profiles is contributed by photons 
with much larger wing scattering numbers. Therefore, a large number of wing scatterings 
acompanied by spatial transfer preferentially in the normal direction may cause the electric 
vectors of the line photons to become perpendicular to the slab normal. From this transfer 
process, we obtain positively polarised fluxes in accordance with our definition of polarisation 
direction. 

Another point we can find in the figure is that the strength of the negative polarisation 
in the far wing regime decreases as ar increases. This effect can be understood in terms of 
the beaming effect. Ahn et al. (2002) investigated the beaming of Lya photons for the cases 
of very large aro, which will also be shown in the next subsection. They showed that most 
photons escape from the extremely thick slab-like media favourably in the direction normal 
to the slab plane, as is the case with an electron scattering thick cloud (Phillips & Meszaros 
1986). If this is the case, the fraction of Lya photons emergent in the direction parallel to the 
slab surface decrease as ar increases. As a result, their contribution to negative polarisation 
diminishes, and eventually the degree of polarisation at the far wings becomes larger and 
eventually positive. Due to the beaming effect, the polarisation flip occurs farther away from 
the line-centre as aro increases. 

The polarisation of Lya in an anisotropically expanding slab was investigated by Lee & 
Ahn (1998). They calculated the cases with the line-centre optical depth r < 10 5 and the 
Voigt parameter a = 4.71 x 10 -4 . Their results show that the polarisation flip occurs near 
the line-centre, and as a result the emergent peaks are negatively polarised on average over 
frequencies. Because the bulk motion of the expanding slab enhances the photon escape, the 
polarization flip occurs at smaller optical depths than in the static medium with the same 
neutral hydrogen column density. 

3.3 Polarisation of the total flux 

Fig. 8 shows the directionality, /(//), which is defined by F(fi)/fi. Here F(fi) is the flux along 
fi. We note that this was also discussed in our previous paper (Ahn et al. 2002). Ahn et al. 
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Figure 8. Directionality of emergent Lya photons for various optical depths to and Voigt parameters a of the scattering media. 
We define the directionality by the flux divided by fj,. The solid thick line represents the limiting behaviour of directionality for 
the Thomson scattering in an opaque electron cloud. We note that the curves is convex upward when aro <C 10 3 in the grazing 
direction. As aro is greater, the directionality gradually converges to the curve for the Thomson scattering in a thick electron 
cloud. 

(2002) discovered that the Lya limb brightening occurs when ar <C 10 3 , and that the Lya 
limb darkening appears when ar > 10 3 . 

In this subsection we show the polarisation (P t ) of total emergent flux. This quantity is 
obtained from the definition P t — J P/(/i)/(/i)/id/i, where -P/(/i) is shown in Fig. 2 and 
is shown in Fig. 8. 

Fig. 9 shows the degree of polarisation as a function of aro- Here aro can be thought to 
be a measure of relative importance between wing scatterings and core scatterings in the 
line transfer. The error in P t is simply propagated from Pf(f^)- The most important thing 
to note in Fig. 9 is that the polarisation sign changes as ar . For ar < 750, the polarisation 
direction is parallel to the slab normal, while for ar > 750 it is perpendiculat to it. 
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Figure 9. Degree of polarisation (Pf) when averaged over both frequency and emergent angle. We show the results for various 
optical depths and Voigt parameters. Pf has the negative value for aro < 10 3 , while the positive values for aro > 10 3 . 

As we saw in the previous subsection 3.1, the case with r = 2 x 10 4 and a = 1.49 x 1CT 2 
has the characteristic optical depth r w = 3.8. We can understand the properties of Lya wing 
scatterings by comparing with the Thomson scattering in an optically thick electron cloud 
with r e t w , because the scattering is described by the same Rayleigh phase function and 
spatial diffusion dominates the transfer process. 

A close examination of the curves in Fig. 2 and Fig. 5 in Phillips & Meszaros (1986) and 
those in Fig. 2 and Fig. 8 in this work shows that when ar < 10 3 the photon flux emergent 
in the grazing direction of the slab surface increases as ar decreases. On the other hand, 
the beaming of emergent Lya is enhanced as ar increases. As a result the critical frequency 
of the polarisation flip moves toward the larger value, which we can see in Fig. 3-7. 
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4 SUMMARY 

We have investigated the linear polarisation of Lja transferred in a static, dustless, and 
thick slab that is uniformly filled with neutral hydrogen. The polarisation behavious of the 
emergent Lja emission can be understood qualitatively by comparing with the polarisation 
behaviour of the Thomson scattered radiation. When the scattering medium is extremely 
thick or ar > 10 3 , photons emerging in the grazing direction are polarised perpendicular to 
the slab normal with the degree up to 12%, which is the same value as that of the Thomson 
scattered radiation in a semi-infinite thick medium (Chandrasekhar 1960). 

The variation of linear polarisation with ar also shows similar behaviour with that of 
the Thomson scattering (Phillips & Meszaros 1986). The linear polarisation of Lja develops 
in the direction parallel to the slab normal when ar < 10 3 , while it becomes perpendicular 
when aro > 10 3 . Our simulated spectropolarimetry of Lja shows the polarisation flip in 
the spectra. The linear polarisation near the line centre is almost zero. The wing parts near 
the line centre are polarised in the direction perpendicular to the slab normal, while the far 
wing parts are polarised in the direction parallel to the slab normal. The zero polarisation 
around the resonance frequency is caused by resonance scatterings just before escape which 
isotropises the electric vector of the Lja photons. The perpendicular polarisation at the far 
wing parts is caused by Lja photons preferentially emerging in the grazing direction. 

The Lja photons constituting the near wing part of the emergent peaks have large wing 
optical depths, and therefore they are beamed in the slab normal direction, which gives rise 
to positive polarisation. As ar increases, the beaming effect prevents Lja photons with far 
wing frequencies from emerging into the grazing direction. Therefore, the location of the 
polarisation flip in the Lya line can be an important indicator of the value of ar . 
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